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Abstract 

We have studied finite-sized single band models with short range pairing interactions be- 
tween electrons in presence of diagonal Fibonacci modulation in one dimension. Two mod- 
els, namely the attractive Hubbard model and the Penson-Kolb model, have been inves- 
tigated at half-filling at zero temperature by solving the Bogoliubov-de Gennes equations 
in real space within a mean field approximation. The competition between "disorder" 
and the pairing interaction leads to a suppression of superconductivity (of usual pairs 
with zero centre-of-mass momenta) in the strong-coupling limit while an enhancement 
of the pairing correlation is observed in the weak-coupling regime for both the models. 
However, the dissimilarity of the pairing mechanisms in these two models brings about 
notable difference in the results. The extent to which the bond ordered wave and the 
^-paired (of pairs with centre-of-mass momenta = tt) phases of the Penson-Kolb model 
are affected by the disorder has also been studied in the present calculation. Some finite 
size effects are also identified. 



I. Introduction 

The competition between electronic correlation and disorder remains one of the prime 
issues of investigation in condensed matter physics during the last few years [1, 2, 3, 4]. 
However, the effect of such disorder is yet to be fully explored in the context of super- 
conductivity. Some earlier experiments showed that in case of some weak coupling super- 
conductors, T c increases with increasing disorder strength, while for some strong coupling 
materials T c is nearly insensitive to the strength of disorder [5, 6, 7]. On the other hand 
very recent experiments observed [8] destruction of superconductivity with increasing 
disorder strength in some low dimensional superconductors. Experiments studying the 
effects of disorder on superconducting A- 15 materials have shown that in these substances 
T c decreases with increasing disorder [9, 10, 11]. Theoretical investigation of properties 
of superconductors in presence of disorder was addressed by Anderson [12] way back in 
1959. He showed that the presence of weak non-magnetic impurities, does not suppress 
superconductivity appreciably. Quite recently [13] the effect of bulk impurity on supercon- 
ductivity has been studied within the negative-?/ Hubbard model with random disorder 
in the site potentials. It showed that superconductivity is suppressed by disorder. More- 
over, it also showed that disorder introduces spatial inhomogeneity in pairing correlation. 
Models like the negative-f/ Hubbard model, which give rise to superconductivity within 
a short range pairing mechanism, have been studied extensively in the context of high T c 
and other exotic superconductors [14]. Therefore, the studies of interplay between disor- 
der and superconductivity within the framework of short range pairing mechanism would 
be of great interest. In this work we describe a study on the effect of Fibonacci- modulated 
disorder on two different models of superconductivity which support short range pairing. 

The models of short ranged pairing that we have focused on in the present work are the 
negative-f/ Hubbard model [14] and the Penson-Kolb (PK) model [15] respectively. These 
two models have been extensively studied [14, 15, 16, 17, 18, 19, 20, 21, 22, 23, 24] in the 
recent past owing to their tentative relevance in the field of high-T c cuprates and organic 
superconductors. However, the pairing mechanism of the two models are of different 
physical origin. The negative-?/ Hubbard model supports a short ranged pairing due to 
an on-site attraction in sharp contrast with the non-local pairing mechanism generated by 
a pair hopping process in the PK model. Such a non-local pair-hopping mechanism gives 
rise to a very rich phase diagram of the PK model [20, 21, 22, 23, 24] as compared to the 



in presence of "disorder" is expected to reveal a qualitative difference in the nature of 
competition between disorder and pairing correlation. 

The present study has been restricted to one dimension (1-d) because of the possibil- 
ity of checking the present results against the earlier studies on the two models, the PK 
model in particular, by a variety of techniques in 1-d. At least in the asymptotic cases, 
some of rigorous results are always available in 1-d rather in higher dimensions. This led 
us to choose a typical 1-d model of diagonal aperiodicity, namely, a Fibonacci modulated 
sequence of the site potentials, for observing the effect of disorder on these superconduct- 
ing models. This type of quasi-crystalline "disorder" not only interpolates between the 
extreme cases of full grown order and random disorder, but also qualifies for the scope of 
experimental investigations owing to the availability of various quasi-crystalline superlat- 
tices in recent times [25]. Thus our main objective, in this paper, is to understand the 
qualitative manner in which the quasi-crystalline disorder modifies the superconducting 
correlation in two specific models with different pairing mechanism and how does the 
nature of this competition change from the weak- to the strong-coupling limit. 

Our investigations for both the models concentrate on decoupling of the Hamiltonian 
within a mean field approximation (MFA) followed by a self consistent solution of the 
Bogoliubov-de Gennes (BdG) equations in real space [13] for the decoupled Hamiltonians. 
The use of a mean-field approach is usually questionable in low dimensions. However, even 
in low dimensional system, such a technique works satisfactorily in a broken-symmetry 
phase [28]. The two models of short-ranged pairing that we have studied here are known 
to exhibit several broken-symmetry phases in the ordered limit. Previous mean-field 
calculations [14, 24] captured these phases satisfactorily in low dimensions and were found 
to compare well with the results obtained by other methods [19, 21, 22]. Moreover, in 
the present scheme of calculation relevant parameters are determined self-consistently for 
each and every site separately which captures successfully the spatial fluctuations induced 
by the Fibonacci disorder [4]. In fact the present scheme of calculation has already been 
tested with success for a randomly disordered Hubbard model in 2-d [13]. 

In section II we define the Hamiltonians and their generic features. In section III we 
furnish the BdG equations. Section IV deals with results that have been obtained in the 
present calculation. Section V summarizes this work. 
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(a) The negative-f/ Hubbard model: 

The negative-?/ (attractive) Hubbard Hamiltonian with Fibonacci modulation in the 
site potentials is given by: 

H-u = Y,( e i - - (* S c iaCi+i<x + h.c.) + U^2 cItwW* (!) 

where, cj^c^) creates (destroys) an electron of spin er(er =|, I) on the i-th site. n icr = 
c\ a Ci a and n» = + n^. e« is the site energy at the i-th site; it takes on the value €a 

or e B according to the Fibonacci sequence: AB AAB AB AAB AAB AB AAB AAB ; /x is 

the chemical potential and t is the single particle hopping integral. The last term is the 
on-site Hubbard interaction. We will take negative values of U in our calculations for the 
attractive Hubbard model. 

The attractive Hubbard model without any modulation in site potential = 0) has 
been extensively studied [14, 16, 17, 18, 26]. In this limit = 0), there is a competition 
between the single-particle hopping (t) and the Hubbard correlation (U < 0). The on- 
site Hubbard correlation favours formation of localized singlet pairs of electrons while the 
hopping term tends to break the pairs. Due to the local pairing mechanism superconduct- 
ing (SC) state and charge density wave (CDW) phases become degenerate in the ground 
state for a half-filled band at zero temperature [19]. This model has been extensively used 
in describing high-T c and other related superconducting systems [14]. Effect of random 
diagonal disorder on the attractive Hubbard interaction has already been studied [13, 27]. 
It is interesting to observe how the Fibonacci modulation alters its superconducting prop- 
erties. 

(b) The Penson-Kolb model 

The PK model with Fibonacci-modulated site potentials is written below: 

U PK = Y^(ti ~ V) n i - (* S C L C i+l<r + h.C.) ~ VQ^ c\ +n cJ +u Q ^ + h.C.) (2) 
i icr i 

The first two terms have similar implications as in the Hubbard Hamiltonian (1). The 
third term is the pair hopping term which is responsible for transfer of a singlet pair 
of electrons (fj,) between neighbouring sites. V is the nearest neighbour pair hopping 
amplitude. PK model with q = 0, favours formation of singlet pairs in real space due 
to this short range pair hopping and therefore, shows a non-local pairing mechanism. In 
this sense the study of this model is complementary to the study of the on-site pairing 



generalizations have been widely studied over the years [20, 21, 22, 23, 24]. The ground 
state phase diagram of this model for a half-filled band in one dimension is understood 
to a great extent as it stands at present. For V > 0, this model shows superconducting 
instability which corresponds to usual pairing with pairs of zero centre of mass momenta. 
However for V < 0, in the strong coupling limit (beyond some value of V, say V c ), this 
model is shown to exhibit 77-pairing with center of mass momentum q = n [21, 23, 24]. 
There is a phase which does not support SC ordering for V c < V < 0. A real space RG 
calculation showed the existence of a CDW phase in this regime [21] and did not consider 
the possibility of the existence of antiferromagnetic bond order wave (BOW) which was 
later shown to coexist with the CDW within a Bosonization calculation [22]. A mean 
field study, however, showed only the antiferro-BOW state in this region [24]. It seems 
rather interesting how these properties are modified in presence of disorder. We focus our 
attention to a Fibonacci modulation in as in the case of the negative-?/ Hubbard model. 

III. The calculation within the MFA and the BdG equations 

(a) For the negative-f/ Hubbard model with Fibonacci modulation 

We solve the negative-?/ Hubbard model first by decoupling the Hamiltonian in favour 
of superconductivity and then by solving the BdG equations of motion in real space in a 
self consistent manner. The decoupled Hamiltonian looks like: 

H-U = XX 6 * ~ ^) n i - ( f C laCi+la + h.C.) + U < n i\ > C ll C H 
i icr i 

+UJ2< nn > c 1t c *T ~ UJ2 < Tin >< nu > +UJ2< c}^ > Qj.q t 

it i 

+uj2< Caen > c\Ai - U H< c \Ai >< > ( 3 ) 

i i 

The BdG equations of motion for the operators and cj^ are: 

*ct T = [c\ v H-u] = te]_ 1T + tc\ +n - UA\ Cil - U < mi > cj t - (ei - /i)cf T (4) 

ion = [c n , H-u] = -tc\ +n - tc\_ n - UA^ + U < n A > c iX + (e ( - /i)c;| (5) 

where c ia = ^jf- and A- =< cLct >. We study the Fourier transform of the supercon- 
ducting gap parameter as defined by, 

A ff = (l/JV) $>*■*■< A, (6) 

icr 

where, r\ is the position of the i-th site. 



The decoupled PK Hamiltonian is given below: 

(f + V r p{ T )5^cJ i c i+u + h.c. 

i 

-v(E A 4+iA+n + h - c -) - V (E + h - c -) + (^Ja j+1 + c.c.) 

i i 

HVpj^ + c.c.) (7) 

where Aj =< cf^ > and =< cj +1(T c ic7 >. 

The BdG equations corresponding to the PK Hamiltonian for the operators Cjj and 
in real space are given below: 

ic% = (t + ^)cJ_ 1T + (t + ^p n )cJ +1T + VAl_ lCil + VA\ +1 c n - (e t - /i)4 T (8) 
ic n = -(t + Vp\y i+n - (t + V Pi ^)4_ n + ^A 1+1 cJ T + + ( Ci - /x)c n (9) 

We solve the BdG equations in a self consistent fashion to determine the Fourier 
transform of the superconducting gap (6) together with the bond order parameter in 
g-space as given by: 

B q = (l/2N)J2°e iqri Pi*- (10) 



(* + ^! i )£ct T c i+ i T + /i.c. 



IV. Results 

(a) For the Fibonacci-modulated negative-f/ Hubbard model 

The negative-f/ Hubbard model has been studied at half-filling on a one dimensional 
chain for system sizes iV=34 and 144. We have chosen = always. The plot of A q 
in Fig.l reveals a maximum at q — 0. This bears a signature of normal superconducting 
phase with singlet pairs having zero centre of mass momentum. The nature of competition 
between the disorder and correlation strongly depends on the value of U. This can be seen 
from the plots of Ao against U (negative values) in Fig.s 2(a) and 2(b) for N = 34 and 
144 respectively. It is observed that for lower values of \U\, value of A is enhanced due to 
Fibonacci modulation. This is due to formation of double occupancies at the sites of lower 
energy. A crossover takes place at intermediate values of \U\. After that A decreases 
with increasing Fibonacci modulation. This happens due to increased backscattering in 
presence of Fibonacci "disorder". Another important fact to note is that the value of 
\U\ where the crossover takes place becomes lower with increasing system size; this is a 



plot of Aj against i at U — —2.8 (a value of \U\ below which the crossover has taken place) 
for iV=34. This plot clearly depicts the spatial inhomogeneity in pairing brought about 
by the Fibonacci modulation. Such results are in agreement with [13]. In the periodic 
limit the distribution of Aj in space is more or less homogeneous except for the boundary 
effects. In presence of disorder Aj closely follows the Fibonacci pattern. Competition 
between disorder and correlation is further revealed in Fig. 4 showing the plot of A as 
a function of e#. It is interesting to observe that for low values of \U\ (U = —1.0), A 
increases with Fibonacci disorder while for strong \U\ (U = —4.0) A sensibly decreases 
with increasing disorder. This is precisely what we have observed in Fig. 2. 
(b) For the Fibonacci-modulated Penson-Kolb model 

The one dimensional PK model also has been studied for half-filling and for the same 
system sizes (iV=34 and 144). Let us first discuss the case of V > 0. Fig.s 5(a) and 
5(b) show the plots of A q for V—1.0. The peaks of A q at q = for V > reveal that 
the system is in a normal superconducting phase. The effect of Fibonacci disorder for 
the high and low values of V can be seen from Fig. 6, the plot of A against e B at 
V = 1.0 and V = 2.0 respectively for iV = 34. For low values of V, disorder makes 
the value of A increase while the opposite phenomenon occurs at larger values of V. 
However, the fall in A with e B for V = 2.0 is extremely slow compared to the notable 
variation of A for V = 1.0. Thus, even a very strong disorder does not seriously affect the 
superconducting properties of the system at high values of V. This is in sharp contrast to 
the case of large \U\. This can be qualitatively understood in the following way. A closer 
look into the eqn.s (4) and (5) reveals that the Hubbard interaction directly renormalizes 
the "effective" site potentials for the electrons of opposite spins to different extents. Thus 
for the larger values of \U\ the site potential seen by an up-spin electron at a particular 
site differs considerably from that seen by a down-spin electron. Therefore, the possibility 
of formation of singlet pairs is reduced due to the competition between strong correlation 
and disorder. But this should not be so drastic in case of the PK model as suggested by 
eqn.s (8) and (9). Let us now discuss the variation of the superconducting gap A as a 
function of the pair hopping amplitude V (for its positive values). Fig.s 7(a) and 7(b) are 
the corresponding plots for iV=34 and 144 respectively at three different values of the site 
energy e B =0, 0.4 and 1.0 at half-filling. In both the plots it is observed that A remains 
vanishingly small until a value of V, say V , is reached. The value of V decreases with 



larger than the gap in the spectrum due to V, for V < V . So the system is unable to 
realize the effect of V in this region and hence A remains almost zero. For N = 144 the 
finite size gap reduces; consequently the value of V goes down. For V > V , A increases 
smoothly and sharply with V. The interesting point to note is that both in Fig.s 7(a) and 
7(b), A has higher values for e B ^ (compared to the case of e = 0) in the regime of 
small V. In this region the effect of pair hopping is rather weak to form a large number 
of singlet pairs. However, the Fibonacci modulation generates sites of lower energy which 
favour formation of "doublons". As a result, a higher value of A is observed for e# ^ 
than for e B = 0. A crossover takes place at a certain V after which the entire process 
reverses. For large values of V, pair hopping process can generate a large number of pairs. 
But the pairing becomes suppressed in the presence of the Fibonacci modulation. This 
is qualitatively similar to the negative-?/ case. However, the degree of suppression of A 
due to disorder in the PK model is not very high. Consistent with this observation we 
also find that the spatial inhomogeneity in Aj is less pronounced in case of the PK model 
(Fig. 8) as compared to the negative-?/ model. 

Next we consider the case V < 0. We study the superconducting properties for 
larger values of |V|. The plots of A q against q show peaks at q = n in Fig.s 9(a) and 
9(b) suggesting that the system is in 7/-paired (with centre of mass momentum q = ir) 
superconducting phase [21]. The peaks are suppressed by Fibonacci aperiodicity. In Fig.s 
10(a) and 10(b) we show the plots of A w against V in the regime V < — 2 for A=34 and 
A=144. It is found that for V = —2.0, A^ has higher value in the periodic limit (e B = 0) 
for both A = 34 and A = 144 than in the cases of e#=0.4 and 1.0. It is interesting to 
note that the suppression of 77-pairing due to disorder becomes less pronounced in larger 
system size. Consequently the decrease of the r/SC correlation across |V| as 2.0 is rather 
gradual (Fig. 10(b)) with increasing disorder in larger systems. 

Next we discuss the intermediate region of — 2 < V < 0. In this regime our result 
shows a bond-order wave (BOW) phase [22, 24] in the periodic limit (for \V\ < 1.9) which 
is evident from the plot of B q against q (Fig. 11). The bond order parameter shows peaks 
at q = 7r for V = -1.0 in Fig.s 11(a) and 11(b) for N = 34 and N = 144 respectively. 
Interestingly enough, the introduction of aperiodicity enhances the peak in case of A=34, 
whereas for A=144, the peak in the periodic limit is a little bit higher than that in the 
Fibonacci-modulated case. This appears to be a finite size effect. The height of the peak 

/• il T-*/~\TT7-\*P 1i • • I T 7" I / 1 1 1 I T 7" I _ . -1 r\ 1 



which ^-pairing takes place). Fig.s 12 and 13 depict simultaneous variations of the BOW 
order parameter and the 77-paired superconducting gap parameter as functions of the pair 
hopping interaction V in the regime — 2 < V < 0. It is observed in Fig.s 12(a) and 13(a) 
that in the periodic limit the transition from the BOW phase to the 7/-paired phase occurs 
at around \V C \ w 1.9 for both iV=34 and 144. This value of \V C \ increases with increasing 
Fibonacci modulation in case of iV=34. But this does not happen for the case iV=144 as 
observed in Fig. 13. This is a result of the fact that Fibonacci sequence can significantly 
disrupt and delay the formation of ^-paired ordering for smaller system sizes, a fact al- 
ready apprehended in Fig.s 10(a) and 10(b). For smaller system sizes the finite size gaps 
are large enough to mask the effects of the pair hopping which reduces the BOW in the 
region — 2 < V < 0. Therefore the quasiperiodic disorder makes the value of \V C \ shift in 
an appreciable manner. For iV=144, the finite size gap becomes much smaller and the 
effect of the pair hopping dominates. Therefore, \V C \ becomes insensitive to the effect of 
the disorder. 

IV. Conclusion 

Summarizing, we have studied the finite-sized negative-f/ Hubbard and Penson-Kolb 
models with Fibonacci-modulated site potentials at half filling in one dimension. In the 
present work, we used a mean field superconducting decoupling followed by self consistent 
solutions of the corresponding BdG equations in real space. We have calculated the 
superconducting gap parameters for normal and 77-paired phases and the BOW order 
parameter for the antiferromagnetic bond-ordered phase in case of the PK model. For the 
negative-f/ Hubbard model we obtain the normal superconducting phase corresponding 
to centre of mass momentum q — 0. Here we observe a crossover at intermediate values 
of \U\. Below this crossover the pairing correlation increases due to Fibonacci disorder 
while above it the pairing is suppressed by disorder. In case of the PK model we observe a 
similar crossover in the positive V sector. However the suppression of SC correlation due 
to "disorder" is much reduced in the PK model as compared to the negative-f/ Hubbard 
model. In case of the negative-?/ Hubbard model the real space pairing amplitude shows a 
spatial inhomogeneity which closely follows the underlying quasiperiodic pattern while the 
feature is less prominent in case of the PK model. The apparent reason behind this (as we 
have discussed earlier) can be understood from the difference in the nature of competition 



mechanism in the PK model is already known to give results which are qualitatively 
different from the negative-f/ Hubbard model [29]. However, in the present study, the 
qualitative difference between these two has been further clarified via the mechanism of 
competition between disorder and pairing term. In the 77-paired phase of the PK model, 
superconductivity is suppressed by disorder in general. In identifying the bond ordered 
phase our calculations match with predictions of previous calculations [22, 24]. It turns 
out that the transition point from BOW to i]-SC phase is not severely affected by the 
presence of disorder in large systems. 

The present study, thus, reveals the qualitative feature of the competition between 
Fibonacci disorder and short ranged pairing interactions within the framework of the 
negative-f/ Hubbard model and the PK model. It is interesting to note how the micro- 
scopic mechanism of pairing modifies the nature of competition in going from the weak- 
to the strong-coupling regime. It is to be noted here that these set of results do not show 
marked variation in going from N = 34 to N = 144. On the other hand, the finite size 
effects that have been identified in the present calculation are essentially controlled by (as 
discussed in Sec. IV) the finite size gap which goes roughly as 1/N. Thus no dramatic 
change in these results are expected to take place beyond N = 144. It may, however, be 
noted that for obtaining very precise quantitative estimate of quantities, e.g. the transi- 
tion point from the BOW to the 77-SC phase in the PK model, a finite size scaling would 
be necessary. Some future work may proceed in this direction. 

For further investigation in this scheme, it would be interesting to include the repulsive 
Hubbard interaction term in the PK model and see the effect of quasiperiodic modulations 
in such systems. Also a detailed study of the effect of band filling away from the half 
filled case would be very much interesting. The possibility of formation of a CDW [21, 22], 
which is left out in the present calculation, could also be investigated. It would also be 
significant to study the case of finite temperatures and higher dimensions. The effect of 
random diagonal disorder in the PK model also needs some attention. 
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Figure Captions 



Fig.l. Plots of A q , the Fourier transform of the superconducting gap parameter, as a 
function of q for different "disorder" strengths (e B ) in the Fibonacci-modulated attractive 
Hubbard model (for U = -2.8) for (a) N = 34 and (b) N = 144. Scale of energy is given 
by t = 1.0. 

Fig. 2. Plots of the superconducting gap parameter A (for pairs with zero centre-of-mass 
momentum) vs. U in the Fibonacci-modulated attractive Hubbard model for (a) iV = 34 
and (b) N = 144 for different "disorder" strengths (e B ). Scale of energy is given by t — 1.0. 

Fig. 3. Plots of the local pairing amplitude A, vs. site index % for the attractive Hub- 
bard model with different strengths of Fibonacci modulation. The plot clearly reveals 
an underlying Fibonacci pattern which is strengthened for larger values of e B . Scale of 
energy is given by t = 1.0. 

Fig. 4. Plots of A vs. e B in the Fibonacci- modulated attractive Hubbard model for 
different values of U. The seemingly different behaviour for smaller and larger values of 
U can be understood from the crossover noted in Fig. 2. Scale of energy is given by t — 1.0. 

Fig. 5. Plots of A q vs. q in the Fibonacci-modulated Penson-Kolb model for the posi- 
tive V vector for (a) N = 34 and (b) N = 144 for different strengths of "disorder" (e B ) 
(scale of energy is given by t = 1.0). 

Fig. 6. Plots of A vs. e B in the Fibonacci-modulated Penson-Kolb model for differ- 
ent values of pair hopping interaction (V) in the positive V sector. Effect of "disorder" 
is found to be distinctly different for low and high values of V. However, the fall of A in 
the strong-coupling limit (e.g. for v = 2.0) is notably slower compared to the case of the 
attractive Hubbard model (Fig.4) (scale of energy is given by t = 1.0). 

Fig. 7. Plots of A vs. V (positive V sector) in the Fibonacci-modulated Penson-Kolb 
model for (a) iV = 34 and (b) iV = 144 for different strengths of disorder (e B ) (scale of 



Fig. 8. Plot of Aj vs. % for the Penson-Kolb model (positive V sector) with different 
strengths of Fibonacci "disorder" . Fibonacci pattern in A; is very much suppressed com- 
pared to the case of the negative-f/ Hubbard model (Fig. 3) even for a reasonably large 
disorder (e B = 1.5) scale of energy is given by t = 1.0). 

Fig. 9. Plot of A q vs. q in the Fibonacci-modulated Penson-Kolb model for the nega- 
tive V sector for (a) iV = 34 and (b) iV = 144 for different strength of "disorder" (e B ) 
(scale of energy is given by t = 1.0). 

Fig. 10. Plots of the superconducting gap A^ (for r) pairing center of mass momentum = 
7r) vs. V in the Fibonacci-modulated Penson-Kolb model (in the negative V sector) for 
(a) N = 34 and (b) N = 144. Effects of different strengths (e B ) of "disorder" are shown 
(scale of energy is given by t — 1.0). 

Fig. 11. Plots of B q) the Fourier transform of BOW order parameter, is shown as a 
function of q for (a) iV = 34 and (b) N = 144 for the Fibonacci-modulated Penson-Kolb 
model (scale of energy is given by t = 1.0). 

Fig. 12. Simultaneous plots of B w and A n vs. V in the negative V sector of the Penson- 
Kolb model for (a) e B = 0.0 and (b) e B = 0.4 for N = 34. Transition from BOW to 
^-paired phase is clearly visible and the transition point shifts towards larger values of 
\V\ for the increase of the "disorder" strength (scale of energy is given by t — 1.0). 

Fig. 13. Simultaneous plots of B w and A n vs. V in the negative V sector of the Penson- 
Kolb model for (a) e B = 0.0 and (b) e B = 0.4 for iV = 144. The transition point 
(corresponding to the transition from BOW to ^-paired phase) is virtually insensitive to 
the "disorder" strength (scale of energy is given by t — 1.0). 
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